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ABSTRACT 

Here, we present numerical simulations of magnetic flux buoyantly rising from a 
granular convection zone into the low corona. We study the complex interaction of 
the magnetic field with the turbulent plasma. The model includes the radiative loss 
terms, non-ideal equations of state, and empirical corona heating. We find that the 
convection plays a crucial role in shaping the morphology and evolution of the emerging 
structure. The emergence of magnetic fields can disrupt the convection pattern as 
the field strength increases, and form an ephemeral region-like structure, while weak 
magnetic flux emerges and quickly becomes concentrated in the intergranular lanes, 
i.e. downflow regions. As the flux rises, a coherent shear pattern in the low corona is 
observed in the simulation. In the photosphere, both magnetic shearing and velocity 
shearing occur at a very sharp polarity inversion line (PIL). In a case of U-loop magnetic 
fleld structure, the field above the surface is highly sheared while below it is relaxed. 

Subject headings: MHD — Sun: interior — Sun: atmosphere 



Introduction 



It is widely accepted that magnetic fields play an important role in driving solar and helio- 
spheric activity. Observational evidence such as Hale's polarity law and Joy's law of tilts suggests 
that active region magnetic fields take the form of J7-shaped loops of magnetic flux anchored to a 
strong, toroidal layer of flux at the tachocline - the thin interface layer between the radiative interior 
and the convection zone, where active region magnetic flelds are believed to be generated and stored. 
Over the past decade, theoretical investiga tions and numerical modeling seem to have validated 
this simple picture (see e.g., the review by iFanI l2009l and referenc es therein). While the physics 



of flux emergence through the deep interior is well-studied (see e.g. iFan et al.lll993l : ICaligari et al 



I995I : lAbbett et al.ll2000l : lFanll2008l ). the physics of how these magnetic structures emerge through 
the upper convection zone, through the visible surface, and into the Sun's atmosphere is not yet 
well understood. 
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One way to make progress is to use magneto- hyd rodynamic (MHD) models to study different 
aspects of the emergence process. IShibata et al.l (|l989l ) used a two-dimensional (2-D) MHD simula- 
tion to study the rising and expansion of an isolated magnetic loop stru cture formed via a magnetic 
buoya ncy i nstability in a two-tempera ture, layered model atmosphere. lEmonet &: Moreno-Insertis 
(11993) and loorch fc Nordlundl (|l998l ) used both 2-D and three-dimensional (3-D) MHD models 
respectively to study the buoyant rise of magnetic flux tubes through a stratified model convection 
zone. In these studies, a random or twisting component of the ma gnetic field along the tub e was 



necessary for the structure to remain coherent during its ascent. lAbbett &: Fisher I (|2003l ) used 
simulations of twisted fi-loops that had risen coherently through the deep layers of the convection 
zone to drive a 3-D MHD model of the solar atmosphere. In this study, the forces resulting from the 
loops rise through the deep interior played an important role in the dynamic evolution of emerging 
flux in the overlying atmosphere, even though the model corona evolved to a relatively force-free 
state. Recent 3-D MHD simulations of flux emergence in a computational domain containing both 



the upper convection zo ne and corona (|Fan 



2001 



Magara fc Longcopell2003l : lArchontis et al. 



2004 



Manchester et al.l 12004') also illustrate the im portance of surface forces during the emergence pro- 
cess. In particular, [Manchester et al.l (|2004l ) found that magnetic tension forces drove shear flows 
and fac ilitated the flux tube's em ergence and partial eruption. This shearing process was discov- 
ered in 



m 



Manch 



Manchester] (|200l|). 



ester &: Low! (|2000l ) and first illustrated with time-dependent numerical simulations 



While much can be learned from these models, they are limited by their simplified way of de- 
scribing the complex energetics of the solar atmosphere. For example, with an artificially-imposed 
thermodynamic stratification or an adiabatic energy equation, the plasma contained in an emerging 
flux rope will cool dramatically as the flux tube expands, and will unphysically inhibit the tube's 
emergence. In light of this, and other limitations inherent to idealized models, it is desirable to 
include ir i the models t he additional physics necessary to improve the realism of the simulations. 
Recently, lAbbettl (|2007l ) implemented a means of approximating optically- thick radiative transfer 
into a 3-D MHD model whose computational domain contained both the upper convection zone 
and corona. He used this model to simulate the quiet Sun magnetic field, ari d studied the physics 



of small-scale flux emergence and submergence. ICheung et al. I (|2007l . |2008| ) solved the radiative 



transfer equation in local thermodynamic equilibrium (LTE) along with the MHD conservation 
equations, and studied flux emergence in a domain containing the upper c onvection zone and pho- 



tosphere. The flux emergence study of iMartinez-Sykora et al.l (|2008l . l2009l ) is particularly notable, 
as they solve the non-gray, non-LTE radiative transfer equation in their 3-D model, and are able 
to extend their computational domain into the corona. Interestingly, they found that as the flux 
tube emerged, chromospheric plasma rose and formed a high-density structure that was supported 
by the expanding magnetic field. 

In each of these studies, it is clear that surface convection driven by radiative cooling dramati- 
cally impacts magnetic structures both above and below the visible surface. We therefore set out to 
improve the treatment of the energy equation in the 3-D MHD model with the Block Adaptive- Tree 
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Solar -wind Roe Upwind Scheme (BATSRUS) developed at the University of Michigan (j Powell et al 



19991 ). In this paper, we describe our improvements to the code, and use this new treatment to 



study the physics of flux emergence in a combined convection zone-to-corona system. We focus on 
the effect of turbulent convection has on the emerging structure, and whether the shear flows de- 
scribed in [Manchester et al.l (|2004l ) continue to be the principle drivers of magnetic flux emergence 
and energy transfer in a radiatively dominated regime. 

The remainder of this paper is organized as follows. The MHD equations and our modifications 
to the code are discussed in section [2j Step-by-step details of our simulations are given in section 
[3l including initial setup, and the parameters for each run. In section [3] we also present the results 
of our simulations and subsequent analysis. Finally, in section U we discuss the implications and 
significance of our results and their relevance to coronal mass ejection (CME) initiation. 



2. Numerical methods 



2.1. MHD equations 



Within BATSRUS, we solve the MHD equations in conservative form using a second order- 
accurate Roe solver on a block-adaptive Cartesian grid using the seven characteristic waves scheme 
(jSokolov et al.ll2008l ). With source terms for radiation, coronal heating, and optional damping of 
vertical flows, the MHD equations take the following form: 
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dt 
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where p, u, e, p, B are the mass density, velocity, total energy density, plasma pressure and magnetic 
field respectively, g is the gravitational acceleration, which is assumed constant, and r represents 
a time scale over which the optional artificial vertical damping is applied. This term is included 
in the system, since it often proves useful to suppress fast moving vertical shocks in the upper 
atmosphere during the relaxation process, significantly speeding up the process (see S ection | 
Qe includes the additional energy source terms, which can be written as (jAbbettl 120071 ): 



Qe Qrad ~l~ Qcr ~l~ Qdamp- 



(5) 
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Here, Qrad is the optically thin radiative loss term; Qcr is an empirical coronal heating term, and 
Qdamp is the energy loss due to velocity damping. In the solar corona, high-temperature low-density 
plasma dominates. Radiative losses in the optically thin limit can be expressed as 



Q 



rad 



-nenpk{T), 



(6) 



where and Up are the electron and proton number densities. Af T) is the total radiative cool- 
ing curve, calculated from Version 6.0 of the CHIANTI database ([Pere et al.l 119971 . |2009| ). We 
approximate optically-thick surface cooling by artificially extending the cooling curve to lower 
temperatures. The temperature and density cutoffs are determined by calibrating our resulting 
stratification against more realistic s imulations o f magneto-convection, where the radiative transfer 
equation in LTE is solved in detail (|Bercikl 12002 ). 



Although we do not have a complete picture of the heating mechanism in the solar corona yet, 
Pevtsov et al.l (|2003l ) demonstrated an empirical relationship between coronal X-ray luminosity and 
unsigned magnetic flux measured at the visible surface. This result, coupled with an assumption 
that the heat in the corona is deposited where the magnetic field is strong, allows us to derive an 
empirically-based coronal heating function of the form 

c(h'^\B\ 



Qc 



jBdV 



(7) 



Here, (/> is the unsigned magnetic flux at the photosphere; B is the magnetic field strength; J BdV 
represents the magnetic field strength \B\ integrated over the region above the model phot osphere; 
and c and a are constants with values of 89.40 and 1.15, respectively, in CGS units (see lAbbett 
20071 ). While we realize that heating due to electron thermal conduction along the magnetic field 
is a crucial contributor to the energy balance of the corona, our focus here is on the dynamic 
emergence of magnetic fiux lower in the atmosphere. We do not attempt to realistically model the 
thermodynamics of the corona, thus we choose to neglect the effects of thermal conduction in order 
to allow for a more thorough exploration of parameter space. This is a limitation of our current 
treatment that will be addressed in future work. 

The third energy source term Qdamp is related to the vertical damping, which can be switched 
off when a relaxed state is achieved in the atmosphere. Since we are solving for the total energy 
equation in the model, the damped energy needs to be taken into account: 

Qdamp — • (8) 



Near the solar surface, 2/3 of th e enthalpy flux is in the fo rm of ionization energy flux, 3 times 
larger than the thermal energy flux (jStein &: Nordlund I l200d ) . Thus, in order to model solar- like 
convective turbulence, we must use a non-ideal, tabu lar equation o f state to close the MHD system. 
For these simulations, we use the OPAL repository (jRogersll200d ) and set the abundance ratios to 
solar values. Specifically, we set X = 0.75, Y = 0.23, and Z = 0.02, where X, Y, and Z refer to the 
mass fractions of hydrogen, helium, and other metals respectively. 
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2.2. Background solar atmosphere 

Active region magnetic fields presumably must travel through the convection zone after their 
formation near the sheared tachocline. It is important to understand the interaction of these 
fields with convective eddies durin g their ascent befo re they reach the photosphere where they are 
observed (see e.g., the review by iFisher et al.l 120001 ) . However, local processes near the surface 
also contribute to the observed evolu tion of the magnetic field. For example, the MHD quiet sun 
simulations by iNordlund et al.l (jl992l ) found that most of the generated magnetic field appeared 
as coherent flux tubes in the vicinity of strong downdrafts. The small-scale convective dynamo 
in granules can maintain a disor dered, locally intense magnetic field, which may contribute to the 
formation of ephemeral regions ([Cattaneo et al.ll2003l ). 



To study the effects of convective motion on flux emergence, we need to generate a realistic 
solar atmosphere with a superadiabatic stratification below the model photosphere. We start 
our simulation by building up a long rectangular domain with dimensions of 0.75 x 0.75 x 15 
Mm^, initialized with uniform temperature and density. The domain is periodic in the horizontal 
directions. At the lower vertical boundary, we modify standard symmetric boundary conditions to 
allow for energy inflow. The top boundary is closed to prevent mass inflows. With these boundary 
conditions, we allow the initially isothermal model atmosphere to evolve self-consistently in response 
to gravitational forces and radiative cooling, and relax to a convectively unstable stratification below 
the surface, coupled to a cold, evacuated upper-atmosphere. The solar-like subsurface stratification 
is achieved by adjusting the density and temperature cutoffs of the extended cooling curve and the 
en ergy iiiput at the base of the mod el convection zone until the resulting atmosphere matches that 
of lBercikI (j^Ooi) and kbbettl koO'^ ) as well as possible. 



We then use this solution as the initial state of a 3-D simulation with a domain of dimensions 
of 10 X 10 X 15 Mm^, which is comparable to the size of an ephemeral region on the sun. Within this 
superadiabatic background atmosphere, convective motion then can be initiated by breaking the 
symmetry with a small energy perturbation applied to the subsurface portion of the computational 
domain. In our simulations, we obtain solar-like convective structures with convective granules of 
upward velocity, up to 6 km/s, surrounded by narrow intergranular lanes of downward velocity up 
to 6 km/s, shown by the gray-scale image in the right panel of Figured! To build a magneto- 
convective state, we introduce a horizontal magnetic field, Bx = 1 G, into the region below the 
surface, while other regions are initially field-free. The magnetic energy is less than 0.001% of 
the kinetic energy of the turbulent motion, which dominates the dynamics in this case. Over 
time, the field is stretched and amplified by convective motions, and flux becomes concentrated in 
intergranular lanes, shown by red lines in the right panel of Figure [H without noticeable deformation 
of the granular convection. Since the focus of this paper is not the study of the convective dynamo, 
we also experiment with stronger initial mean fields to more quickly build up semi-realistic models 
of magneto-convection. 



We then put a 1 G vertical magnetic field into the domain to heat the corona by activating 
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Fig. 1. — Left: The average vertical density (solid) and temperature (dashed) stratification. The 
dash-dotted line indicate the height where the density and temperature values (1 x 10~^ g/cm^, 
5730 K) are comparable with photospheric values. Right: The vertical velocity structure at the 
photosphere {z = -2.5 Mm). Red lines show regions with \Bz\ greater than 3G. 



the empirical coronal heating term Qcr, and open the upper boundary to allow waves to travel 
outward. The vertical variation of the relaxed atmosphere is described in the left panel of Figure 
[U which shows average density and temperature plotted as functions of height. At z = —2.5 Mm, 
the density and temperature values (1 x 10~^ g/cm^, 5730 K) are comparable with the observed 
values of the photosphere, and thus, this layer is defined in our simulations as the solar surface. 

After the generation of a realistic convective state and overlying atmosph ere, we ins ert a 
thin, horizontal twisted flux rope into the model's convection zone. Following iFanI (120011 ) and 



Manchester et al.l (|2004l ). we describe the initial flux rope by 

Boe'^'/'''x + qrBoe-'-'/'''0, 



Br 



(9) 



where we set a = 0.3 Mm. q is the twisting factor deflned as the angular rate of held line rotation 
per unit length in the axial direction. The density is depleted in the central section of the flux rope, 
which produces an il-loop structure by buoyancy using: 



(10) 



where A = 1.5 Mm, and the ratio i] is deflned as 

2 -Doe ' 
rj = -- 
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(11) 



which maintains a force balance. 
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To make the total energy unchanged, we correct the pressure by: 

P = Po(l -V)- 



(12) 



3. Results 

Results from our model are presented as follows. In Section I3.H we describe the steps we took 
to simulate the rise of a buoyant flux rope, and the structures that developed during the course of 
its ascent. In Section [3.21 we discuss the shearing motion during the emergence process, and finally 
in Section 13.31 analyzes a particular case of "U-loop" formation. 




Fig. 2. — The 3-D structure of magnetic field of Run 1 (a), Run 2 (b) and Run 3 (c) at t = 20, 
20, 16 min, respectively. The planes show the Bz on photosphere. Black lines outline regions with 
downflow speed greater than Ikm/s. Flux ropes with more twist or strength have less horizontal 
expansion and more coherence during emergence. 
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3.1. Emergence of the Flux Rope 

In the sections below, we experiment with different twisting and magnetic strength values 
shown in Tabled] and evaluate their role in affecting the emergence. At the photosphere, the flux 
rope first emerges as two closely-aligned flux concentrations with opposite polarities. Then these two 
polarities start to separate from each other since the buoyant section continues to lengthen. Figure 
[2] shows the 3-D structure of magnetic fields for each run, with the color showing the photospheric 
Bz field, and black lines marking downdrafts with speed greater than 1 km/s. For Run 1, with a 
weak field Bq = 7.0 kG and low twist q = 1.0, we observe little coherency in the emerging flux 
rope, which expands preferentially in the horizontal direction. The horizontal expansion occurs 0.5 
Mm above the photosphere, at the temperature minimum. And on the photosphere, convective 
motion dominates over the rising motion of magnetic flux, as seen in Panel (a) of Figure [2j The 
twisting factor was increased in Run 2, and a coherent flux rope was observed in the lower corona 
after emergence, shown in Panel (b). In Run 3, both twisting factor and flux rope strength are 
increased compared with Run 1. We find that the increase of a non-axial component of the initial 
magnetic field (Panel (c) of Figure [2]) enhances coherency, and prevents the emerged rope from over- 
expansion in the horizontal direction, as shown in Panel (a). Turbulent plasma fiows distort the fiux 
rope during its rise, while a fiux rope with strong mag netic field deforms the convective granules 



after emergence, consistent with ICheung et al. I (|2007l ). However, regardless of the deformation 



the strong vertical magnetic fields concentrate in the downdrafts after fully emerging above the 
surface. In the following sections, we will analyze Run 3 in detail. 



3.2. Shearing Motion 



Photospheric shear ing motions have long been observed to be in s t rong association with solar 
fiares and CMEs (e.g., iMeunier fc KosovichevI bood : Ivang et all bood : ISchriiver et allboosl ). Be- 
sides the photospheric motion, lAthay et al.l (119851 ) found a shearing area in th e transition region 
and upper chromosphere. Helioseismic flow maps (jKosovichev Duvalll l2006l ) also indicate the 
magnetic energy buildup prior to CMEs is driven by the strong shearing flows with speed 1-2 km/s 
at the depth of 4-6 Mm below the surface. 



Table 1. Parameters for runs 



runs 


q 


Bo (kG) 
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-1.0 


7.0 
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Fig. 3. — The evolution of sheared (sohd) and emerged (dashed) energy flux at the photosphere for 
Run 2 (left) and Run 3 (right). In the right panel, dotted and dash-dotted lines show the evolution 
of shear and emerged fluxes, respectively, at z = -0.45 Mm. 



In our simulations, coherent shearing flows develop below the photosphere and extend into the 
lower corona in the later phase of emergence. We calculate the components of the Poynting flux 
associated with horizontal and vertical motions on the photosphere by: 



Eshear = " / — (BxUx + ByUy) B^dS, (13) 



1 

8^ 

Eemerge = f ^ {^1 + B^) U,dS. (14) 



The temporal evolution of the energy fluxes is shown in Figure O In Runs 2 and 3, the top 
of the flux rope reaches the photosphere at times t = 9.5 and 5.5 min, respectively, shown in 
Figure El When the magnetic field first emerges at the surface, a simultaneous increase in sheared 
and emerged energy flux occurs, with a greater increase in emerged flux. The emerging motion 




-4 -2 2 4 -4 -2 2 4 

X (Mm) X (Mm) 



Fig. 4. — The vertical magnetic field structure on photosphere at times t = 10 (a) and 15 (b) min 
for Run 3. Blue lines show regions with |(V x u)^! greater than 1.5 times of the average value. The 
horizontal magnetic field is shown by the red arrows. The locations of strong PILs show a good 
coincidence with regions with higher |(V x \i)z\ values. 
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dominates the energy flux for the first 10 minutes, during which the buoyant part of the flux rope 
evolves to a fuUy emerged state. Afterwards, near the surface, the magnetic pressure gradient force 
is balanced by gravity and gas pressure, which slows down the emergence and decreases the emerged 
energy flux after the saturation of emergence. The emerged energy flux drops to negative values 
in the later phase due to the concentration and subduction of the magnetic field in downdrafts. 
Dotted and dash-dotted lines in right panel of Figure [3] show the shear and emerged energy fluxes, 
respectively, in lower corona, 2.05 Mm above the photosphere, which has a temporal delay of 6 min 
but represents a similar trend as the fluxes at the photosphere. 
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Fig. 5. — The vertical magnetic field structure on photosphere at times t = 10 (a), 12 (b), 15 (c) 
and 20 (d) min for Run 3. The red arrows show the horizontal velocity field. In panel (d), the green 
line outlines the region with velocity shearing illustrated by |(V x u)^| greater than the average 
value. In panel (c), the green box marks the subregion shown in Figure [7] and Figure [H High 
velocity shearing occurs along the PILs during the emergence. 



When the rising motion succumbs to photospheric downfiows, the main energy release into 
the overlying atmosphere is contributed by the shear flow, which continues long after the flux has 
emerged through the photosphere. The sheared energy fiux incre ases by a factor of 4 froin the initial 
emerging phase. The energy evolution resembles that shown in [Manchester et al.l (|2004l ). with the 
fundamental difference that active convection ultimately subducts magnetic field and reverses the 
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emerging energy flux. 
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Fig. 6. — Variation of u^; on X = plane at times t = 5 (a), 10 (b), 15 (c) and 20 (d) min for Run 
3. Black lines show the magnetic field lines on Y-Z plane, with arrows indicating the direction. 
Dashed lines indicate the height of the photosphere. Coherent shearing pattern develops in the 
lower corona and extends down to the chromosphere when the rope rises. 

Figure H] shows the evolution of vertical and horizontal magnetic field structures at the photo- 
sphere for Run 3. The horizontal gradients in velocity are illustrated with contour lines of (V x u)^. 
At t = 10 min, Panel (a) of Figure H] shows a highly sheared magnetic field structure at this initial 
phase. Panel (a) of Figure [5] shows a horizontal velocity field flowing outward from the center of 
the granule, at t = 10 min, when the top of the flux rope reaches the photosphere and the size of 
the convective granule lying above it increases. After that, the emerging motion is saturated, and 
the fragmentations and distortions of the magnetic field mainly result from shearing and convective 
motions. In Panel (b) of Figure \5\ the two polarities start to move along the PIL in opposite 
directions, and form strong PILs shown in Panel (c) and (d). At time t = 15 min, Panel (b) of 
Figure m shows a good coincidence of the PIL and regions with (V x u)^ greater than 1.5 times of 
the average value. This implies a high velocity shearing along the PIL at the photosphere in the 
later phase of emergence, which is clearly seen in the region outlined by the green box in Panel (d) 
of Figure [5l 
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Figure [6] shows the evolution of the horizontal velocity, u^, and the magnetic field lines on the 
X = plane for Run 3. A coherent shearing pattern with velocity Ux up to 20 km/s develops in 
the lower corona and extends down to the chromosphere and p hotosphere during the rising of the 
rope. The shear pattern is similar to [Manchester et al.l (|2004l ) with notable exceptions of being 
more highly structured with reduced velocity. We also find that shear flows in the corona persist 
even when the convective motion washes out any coherent shearing at the photosphere. 
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Fig. 7. — (a) Fy structure at time t = 16 min (b) Fy structure at time t = 18 min (c) y component 
of velocity perpendicular to magnetic fields at time t = 16 min (d) y component of velocity 
perpendicular to magnetic fields at time t = 18 min in region marked by green box in Panel (c) of 
Figure [5] for Run 3. Black lines show the magnetic field lines on the plane, with arrows indicating 
the direction. Dashed lines show the height of the photosphere. The field lines form a U-loop 
structure with a shearing velocity up to it 1.5 km/s on each side of the PIL (outlined by the green 
boxes). 
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Fig. 8. — (a) Ratio of thermal pressure to magnetic energy logiQ{p/ Emag) at time t = 16 min (b) 
ratio of kinetic energy to magnetic energy logio{Ek/Emag) at time t = 16 min. Black lines show the 
magnetic field lines on the plane, with arrows indicating the direction. Dashed lines show the height 
of the photosphere. Kinetic and thermal energy dominate in the domain below the photosphere, 
while above it magnetic energy dominates. 




X(km) 

Fig. 9. — Bottom view of the magnetic field lines in region marked in green box in Panel (c) of 
Figure [5] for Run 3. Plane shows the Bz structure on photosphere. The magnetic field lines run 
perpendicularly to the PIL, minimizing the energy. 
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3.3. U-loop Structure 

Here, we examine a U-loop - a mag netic feature formed when magnetic flux is pulled down 



into an intergranular lane, as reported by iTortosa-Andreu &: Moreno-Insertia (j2009l ). We evaluate 



the role of Lorentz force F in a subregion outlined in Panel (c) of Figure [5] where the U-loop forms, 
by examining Fy and the y component of velocity perpendicular to the magnetic field {u±_)y. This 
velocity component, {u±)y, excludes siphon flows parallel to the magnetic fleld and thus illustrates 
the motion of magnetic field lines. In Figure [71 we show that the field lines on the X — Z plane 
evolve to a U-loop structure with a shearing velocity {u±)y up to ± 1.5 km/s on each side of the 
PIL (outlined by the green boxes). The upper two panels illustrate the magnitude of the Lorentz 
force while the bottom panels show the magnitude of {u±)y. The fact that Fy and {u±)y are in the 
same direction-reversing across the neutral line (outlined by the green boxes) suggests the velocity 
shearing is due to the Lorentz force. Panels (c) and (d) of Figure [5] show a highly sheared magnetic 
field and velocity field in this region on the photosphere. However, when viewed from the bottom, 
as shown in Figure [9] (with 3-D magnetic lines over a gray-scale image of the photospheric field 
in this region), the magnetic field lines are unsheared below the surface, running perpendicularly 
to the PIL, and minimizing the energy. Figure [8] shows the ratio of the thermal pressure and 
kinetic energy to the magnetic energy in the region in Figure [71 Below the photosphere, kinetic 
and thermal energy dominate, while above the photosphere, convective motion stops and gives way 
to the magnetic energy. 

It has been commonly observed that the central regions of f2-loops become highly sheared as 



they emerge and unshear as they submerge (jManchesteii |2007| ) . We find for the first time that 



magnetic U-loops will unshear by fiows driven by the Lorentz force, as illustrated by the gradient 
in By. These fiows transport energy and axial flux from the convection zone to the corona, where 
the U-loop is expanding. This simulation further illustrates the dynamic coupli ng of the convectio n 



zone and corona, due to shearing flows driven by Lorentz force, suggested by iManchested (|2007l ). 
This coupling occurs for two reasons: the transport of axial flux and the requirement that axial 
fleld strength be equilibrated along field lines to achieve a state of force balance. 



Discussion and Conclusions 



Observations have shown that CMEs originate from polarity inversion lines on the photosphere, 
where strong rnagnetic and velocity s hearing occurs. Study of vector magnetograms of active regions 
(|Falconeiil200ll : [Falconer et al.ll2006l ) found the nonpotentiality of the magnetic fields, i.e. the total 
free magnetic energy, to be the determinant of the CME productivity. Other observations also 
indicate that solar eruptive events, e.g. flares, filaments, and prominences, preferentially occur 
in magnetic structures where the n iagnetic field runs parallel to the magnetic inversion line, i.e . 
highly sheared m agnetic stru c tures (Foukaj[l97l[; [Leroy|[l989l : [Canfield et al.|[l999l : [Liu et al.[[2005l : 
Su et al.[ [2007[ ). [Peng et al.[ (|2006[ ) observed persistent and long-lived (> 5 hrs) shear fiows in 
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Active Region NOAA 10486 along the neutral line prior to the occurrence of the XIO flare on 2003 
October 29, which may be explained by the successive emergence of a much larger and stronger 
flux rope/system than that in our simulations. 

Surface plasma flows have long been recognized to have an irnportan t infl uence on the evolu- 



tion o f coronal magnetic fields. Many simulations, e.g. iMikic et al.l (119881 ) and Ivan der Hoist et al 



(|2007l ). have been performed with flux ropes and arcad es whose footpoints are subjected to surface 
motions imposed as boundary conditions. For example. I Amari et al.l (|2003l ) has shown that converg- 
ing motions toward the magnetic inversion line can drive a catastrophic nonequili brium transition 



into e ruptions when accompanied by shearing motions along the inversion line. In lAntiochos et al 



(|l999l ). sheared arcades reconnect with a surrounding flux system, resulting in an eruption, in which 
the magnetic free energy stored in the closed shear arcade is released. In these coronal models, 
shear flows are only pre scribed as bounda r y conditions to trigger i nstabilities and eruptions. Recent 
numerical simulations (jManchesteii l2003l : [Manchester et al.l |200J) have shown shear flows develop 
self-consistently during the emergence of the flux system (arcade and rope respectively) through 
a stratified atmosphere. The shearing motions result from the Lorentz force that occurs as the 
magnetic field expands into the stratified ambient atmosphere. This shearing mechanism explains 
the following: (1) coincidence of magnetic and velocity inversion lines; (2) the time evolution and 
magnitude of the shear flows in different layers of the atmosphere; (3) the large scale shear pat- 
tern, which is most concentrated at the PIL; (4) transport of flux and energy from the convection 
zone to the corona; (5) loss of equilibrium when the shear flows are unable to equilibrate the axial 
flux along magnetic field lines - a process that c an explain many obseryed occurances of CMEs , 



eruptive filaments, and f lares (jManchester 



Archontis fc Torokll2008l : iMacTaggart fc Hood 



2003 



200' 



Manchester et al.l 120041 : iManchested l2007l . l2008l : 

^). However, in these idealized models without 



the convective motion, shear flows grow rapidly in the absence of competing mechanism. In ad- 
dition, the use of a purely adiabatic energy equation leads to expanding flux ropes whose plasma 
cools far too rapidly, thus inhibiting the emergence process. 

In light of these problems, we incorporated additional physics into our MHD models, and 
simulated the buoyant rise of a twisted flux rope as it emerged through a turbulent model convection 
zone into a model corona. Flux ropes with more twist or greater magnetic strength can rise through 
the convection zone coherently, while less twisted or weaker flux ropes tend to expand in the 
horizontal direction, while being distorted by the turbulent motions, as shown by Figure |2j During 
its emergence, shearing motions at the photosphere develop and extend into the lower corona, with 
a speed up to 20 km/s, which is comparable with the shearing velocity of 20 km/s reported by 
Chae et al.l (|2000l ) in active region loops. In our simulations, shearing motion couples the emerged 
flux with the subsurface magnetic fields, and transfers the energy from below the surface into the 
solar atmosphere. In Run 3, where the twisting factor q = —1.5 and magnetic strength Bq = 14.0 
G, the total sheared energy can go up to 1.08 x 10^^ ergs in 21 min, about 1% of t he average amount 



released in CMEs. The consistency between our results and previous simulations (jManchester et al 



20041 ) shows that the physical mechanism of Lorentz force-driven shear flows is very robust during 
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the flux emergence in a domain comparable to ephemeral regions, regardless of the presence of 
thermodynamic processes such as turbulent convection. Given that active regions are 10 times 
(with 100 times the area) larger than this simulated flux concentration, our results suggest that 
shear flows driven by the Lorentz force can readily provide the free magnetic energy for CMEs. 

In order to study the triggering mechanism for the explosive activities associated with magnetic 
field structures, we will extend our domain size to be comparable with the size of active regions 
and increase the flux rope size. Our model is very useful in the sense that it can be expanded 
in horizontal and vertical directions using the adaptive grids. The model includes different layers 
and physical processes in the solar atmosphere and new features can be easily implemented. In 
future simulations, we need to take into account the background magnetic field and the field- 
aligned heat conduction, and investigate the effec ts of surface footpoint motions on reconnectio n 



between emerging a r id arn bient field, as reported bvlArchoiitis et al.l (|2004l ) , iGalsgaard et al.l ()2007l ) 
Archontis Torokl (|2008l ) and iMacTaggart k HoodI (|2009l ). 
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